Spin Dynamics in Pyrochlore Heisenberg Antiferromagnets 



o 
o 

(N 
< 



I 

O 

o 



o 

o\ 
o 



X 



p. H. ConlorB and J. T. Chalker 
Theoretical Physics, Oxford University, 1, Keble Road, Oxford 0X1 3NP, United Kingdom 

(Dated: April 9, 2009) 

We study the low temperature dynamics of the classical Heisenberg antiferromagnet with nearest 
neighbour interactions on the pyrochlore lattice. We present extensive results for the wavevector 
and frequency dependence of the dynamical structure factor, obtained from simulations of the pre- 
cessional dynamics. We also construct a solvable stochastic model for dynamics with conserved 
magnetisation, which accurately reproduces most features of the precessional results. Spin correla- 
tions relax at a rate independent of wavevector and proportional to temperature. 
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Geometrical frustration in magnets inhibits ordering. 
Simple, classical models for these systems have very de- 
generate ground states [H, 0]. Reflecting this degener- 
acy, highly frustrated magnetic materials characteristi- 
cally remain in the paramagnetic phase even at tempera- 
tures low compared to the scale set by exchange interac- 
tions. Behaviour in this cooperative paramagnetic regime 
has been the focus of much recent research 01 • 

Nearest neighbour antiferromagnets on the pyrochlore 
lattice with classical n-component spins are representa- 
tive of a large class of models [1, 0] . They have remark- 
able correlations at low temperature, which are interme- 
diate between those of conventionally ordered and com- 
pletely disordered systems. These can be understood 
by mapping spin states onto configurations of a vector 
field, or flux field, which is solenoidal for ground states 
0, B 0] ■ Gaussian fiuctuations of this flux field provide 
a coarse-grained description of the cooperative param- 
agnet. Static spin correlations have a power-law depen- 
dence on separation, inside a correlation length ^ that 
diverges as temperature T approaches zero. These cor- 
relations result in sharp features, termed pinch points, in 
diffuse scattering as a function of wavevector. 

The dynamics of cooperative paramagnets has not 
been studied as extensively as the statics, but some in- 
gredients are clear. In a Heisenberg model with preces- 
sional dynamics, the short-time behaviour can be viewed 
in terms of harmonic spinwave fiuctuations in the vicin- 
ity of a specific ground state, while over longer times 
the system wanders around the ground state manifold. 
This second component to the motion results in decay 
of the spin autocorrelation function at long times, with 
a decay rate shown to be linear in T using simulations 
and phenomenological arguments Q. These theoreti- 
cal ideas are supported by inelastic neutron scattering 
measurements on pyrochlore antiferromagnets: an early 
study of CsNiCrFg revealed strong temperature depen- 
dence to the width in energy of quasielastic scattering for 
T < |0cw| while recent work on Y2RU2O7 shows a 



precessional dynamics than has been available so far. The 
topic is interesting from several perspectives. First, in 
view of the pinch points in static correlations, it is natural 
to ask about the wavevector dependence of the dynami- 
cal structure factor, accessible in single-crystal measure- 
ments. Little is currently known about this: the auto- 
correlation function of Ref [H] is expressed as an integral 
over all wavevectors, while the measurements of Ref. [IJJ 
used a powder sample. Second, dynamics in the param- 
agnetic phase of unfrustrated antiferromagnets is domi- 
nated by spin diffusion 13, [13, [13, and one would like 
to know whether this extends to the cooperative param- 
agnct. Third, behaviour in the Heisenberg model should 
be compared to that in spin ice, which is repesented by 
the Isiirg pyrochlore antiferromagnet with dynamics con- 
trolled by the motion of monopole excitations • 

In outline, our results are as follows. We find at 
low temperature three types of behaviour in different re- 
gions of reciprocal space, (i) Close to reciprocal lattice 
points, correlations are dominated by spin diffusion with 
a temperature-independent diffusion constant, (ii) At a 
generic wavevector [not included in (i) or (iii)] correla- 
tions are Lorentzian in frequency with a width linear in 
T and independent of wavevector. (iii) Close to nodal 
lines in reciprocal space on which the static, ground- 
state structure factor vanishes [8], dynamical correla- 
tions are dominated by finite-frequency spinwave contri- 
butions. This picture is hence very different from that 
for the kagome Heisenberg antiferromagnet, which shows 
order-by-disorder and propagating modes [16]. 

We consider the classical Heisenberg antiferromagnet 
with nearest neighbour interactions on the pyrochlore 
lattice. Lattice sites (labelled form corner-sharing 
tetrahedra (labelled a,/3). Spins Si are unit vectors and 
Lq, = X)iea ^4 total spin of tetrahedron a. The 

Hamiltonian is 



c. 



(1) 



ill 



width linear in T as predicted 

Our aim in this paper is to establish a much more com- 
prehensive description of cooperative paramagnets with 



where c is a constant. Ground states satisfy Lq = for 
all a. The equation of motion, describing precession of 
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each spin around its local exchange field, is 



(2) 



where sites j are the nearest neighbours of i. The global 
spin rotation symmetry of Eq. ([1]) implies conservation 
of total spin in the dynamics. 

Before presenting results from a molecular dynamics 
study of Eq. we consider an analytically tractable 
stochastic model for the dynamic behavior. It is known 
that static spin correlators for the classical Heisenberg 
model are well described by those for 7i-component spins 
in the large n limit [1, 0]- Building on this, we set out 
to endow the n = oo model with appropriate dynamics. 
First we recall some details of the static model. Taking 
the second form of the Hamiltonian in Eq. ([1]), a single 
spin component in the large n limit has the unnormalised 
probability distribution e~^^ with 
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(3) 



where now 1^ = X^jgq sum of 'soft' spins 

(— oo < Si < oo) on tetrahedron a. The spin length 
is constrained by the Lagrange multiplier A. For /3 — > oo, 
the second term in Eq. ^ enforces all the la to be zero. 
The interaction term written directly in terms of the 
spins is ^PJJ2iji^ij + '^^ij)siSj where Aij is the ad- 
jacency matrix for the pyrochlore lattice. We call the 
combination Aij + 2Sij the interaction matrix. Its eigen- 
values f^(q) are labeled by wavevector q and a band 
index /i £ {1,2,3,4}. Two bands are flat (wi_2(q) = 0) 



and two (/i — 3,4) are dispersive. Requiring 
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to mimic behaviour of a single spin component in the 
Heisenberg model, A = 3/2 + 0{T/J) for T < J. We 
denote the Fourier transform of the spin variables Si by 
Sq, where a is a sublattice index, and define the sublat- 

tice sum Sq = X]t=i ■^q- Transforming from to the 
basis (denoted by tildes) that diagonalises the interac- 
tion matrix gives collective spin variables s^. We want 
to introduce time dependence and calculate the dynamic 
correlation function S{q,t) = (•Sq(t)s* (0)), and its time 
Fourier transform, the dynamic structure factor S'(q, cj), 
measured using neutron scattering. 

There are many choices of dynamics which repro- 
duce any given equilibrium distribution. To approximate 
Eq. (121) we demand a local dynamics that conserves the 
total spin. We can ensure this by requiring the spin on 
each site to satisfy a local continuity equation. We in- 
troduce spin currents on bonds of the pyrochlore lattice, 
which have drift and noise terms. We take the drift cur- 
rent on a bond linking two site to be proportional to the 
difference in the generalized forces dE/dsi at the sites. 
This favours relaxation towards a configuration that min- 
imizes E] the thermal ensemble is maintained by noise 



which has an independent Gaussian distribution on each 
bond. These assumptions lead to the dynamical equa- 
tions for the soft spins 



dsi ^ A dE ^ , . 



(4) 



where the matrix A is the lattice laplacian (for a lat- 
tice with coordination number z, A^; = An — zdn). 
The correlator of the noise Q{t) at site i, {Ci{t)Cj{t')) = 
2TTAijS{t — t'), has an amplitude fixed by the require- 
ment of thermal equilibrium. The only free parameter 
in the model is the rate F, which sets a timescale for 
dynamical processes. The Langevin equation Eq. (j4]) is 
straightforward to solve in the diagonal basis. It gives 
the correlation function 
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The dynamic correlation function is then 



5(q,t) = ^5M(q)(5(^WS^q(0)> 



(5) 



(6) 



where the structure factors gf^iq) are formed from the 
eigenvectors of the interaction matrix. They satisfy the 
sum rule J2l=i9t^i'i) = 4- 

For completeness we present their explicit forms here. 
In the notation of 0], where Cab = cos ( "^"f and 



with Q 



-3 and 



defining P = ^/l + Q, the eigenvalues of the interaction 
matrix are ui,2 = and V3^4 = 4 =p 2P. Further defining 
si = sin^ (^) and C(^ab) = Cab + c^, the g^{q) are for the 
degenerate flat bands 



g = 91+92 
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3-Q 

and for the dispersive bands 

53,4 - 2-^5 (1 ± 2P-i)±p-i(2-C(j,,)-C(,,)-C(,^)) 

which indeed satisfy 5i -I- 52 + 53 + 54 = 4. 

We now examine the implications of this model for 
T <ti J, emphasizing the features (i) - (iii) mentioned in 
our introduction. From the exponent in Eq. ([5|) we ob- 
tain a characteristic time for decay of correlations, (i) In 
the vicinity of q = only the coefficient 54 (q) is non- 
zero and so behaviour is controlled by the fourth band 
whose decay rate is t^^ = STJa'^q'^ + 0(g**), where a is 
the pyrochlore site spacing; from this we identify the spin 
diffusion constant in this model as D = STJa^, indepen- 
dent of T. (ii) At a generic wavevector where gi and 
(72 are non-zero most of the spectral weight is in the flat 
bands, with decay rate = SFAT, independent of q; in 
an approximation where only the flat bands contribute, 
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FIG. 1: (Colour online) Evidence for spin diffusion, (a) 
Dependence of S{q, u) on q and uj at small q. (b) Scaling 
collapse following Eq. (O at multiple temperatures (/3J = 
20, 40, 60, 80) and four of the wavevectors plotted in (a). Also 
plotted is the prediction of Eq. (|6]). 

this implies the dynamic structure factor factorizes as 
S{q,u>) = 5'(q)/(w), a possibihty noted in (iii) On 
nodal lines |8|, high symmetry directions in reciprocal 
space along which 51 (q) + 52(q) = 0, the decay rate is 
wavevector-dependent and 0{J) away from q = 0. 

To test these ideas we have performed simulations of 
the full precessional dynamics, Eq. Low T config- 
urations are obtained using a Metropolis Monte Carlo 
sampling method. We take these as initial configurations 
for numerical integration of the equations of motion us- 
ing a 4th order Runge-Kutta algorithm with adaptive 
step size. Energy and total spin are conserved to rela- 
tive errors no greater than 10~®. We report data from 
simulations on system sizes with total number of sites 
N = AL^ for L = 16,32. We calculate the dynamic 
correlation function S{q,t) = (Sq(i) • S_q(0)) and the 
dynamical structure factor S'(q, w) = (|Sq(aj)p) where 
Sq = J2t=i is Fourier transformed spin configu- 
ration and a runs over the four sublattices. We present 
results under the headings (i) - (iii) as above, expressing 
q in reciprocal lattice units as q = 2nk for Figs. [H [31 SI 

(i) Since the total magnetization is conserved, one ex- 
pects diffusion at sufficiently small q. The simulations 
confirm diffusive behavior with a diffusion constant in- 
dependent of temperature. At small the data should 
collapse onto the scaling form appropriate for diffusion, 

2D 

/3g^'5'(q,cj) = 3x-^-v^-2— ^ , (7) 

where x is the susceptibility per primitive unit cell and 
the factor of 3 is due to the 3 spin components. Fig. [T] 
shows this scaling collapse when plotted as in Eq. ((7|) , and 
demonstrates that the diffusion constant is independent 
of temperature. The prediction of Eq. whose small 
q limit has precisely the form of Eq. ([7]), is an excellent 
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FIG. 2: (Colour online) Main panel: decay rate of S(q, t) 
as a function of T at wavevectors as indicated. Upper inset: 
S(q,ij) at 13 J = 40 (red) and 100 (blue), for k = (1,1,1). 
Lower inset: ^(q, lo) on a nodal line at k = (0, 0, 1.25). 

fit with F = 0.167. 

(ii) At a generic point in reciprocal space (not near q ~ 
or a nodal line) the structure factor is well described by 
a Lorentzian centered on w = 0, indicating relaxational 
dynamics (Fig. [2l upper inset). The decay rate for this 
relaxation (Fig. [21 main panel) is proportional to T and 
independent of wavevector, even close to the pinch points. 

(iii) On nodal lines, by contrast, the width in frequency 
of S{(i,ijj) is 0{J) and depends little on T: see Fig. [21 
lower inset. High frequency (w = 2. 5 J) and zero fre- 
quency behaviour is also presented in Fig. [31 as a survey 
of S{(\,uj) in the {qx,qx,qz) plane, using data taken at 
/3J — 500 and typical of all low temperatures. Weight 
in S'(q, Lu) at LO J and low T can be viewed as due to 
spinwave fluctuations in the vicinity of an instantaneous 
ground state. In contrast to behaviour in the kagome 
antiferromagnet [l6j , there in no evidence in Fig. [31 for 
sharp propagating modes. 

We next demonstrate that the stochastic model accu- 
rately reproduces most aspects of the precssional dynam- 
ics. We examine behaviour with each type of dynamics 
at different temperatures and times along a path in the 
BriUouin zone: V = (0,0,0) ^ (2,2,2) S (0,0,2). This 
consists of the section Vi along high symmetry nodal 
lines passing through a pinch-point, and a section V2 
typical of reciprocal space. We compare in Fig. [H the dy- 
namic correlation function, normalized to unity at t = 0, 
at two temperatures with the predictions from Eq. ^ 
along the path V . The excellent agreement of the curves 
across multiple temperatures, wavevectors and times is 
good evidence that the stochastic model is sufficient to 
capture the relaxation behavior. It fails only at short 
times {t < J~^, not shown in Fig. [J]) on nodal lines, 
where it does not account for oscillatory spinwave con- 
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FIG. 3; (Colour online) Intensity map of l3S{q^, to) in the 
{qx,<lx,qz) plane at /3J = 500. (a) Lower panel: zero fre- 
quency (data divided by 1.25 x 10^); white line is path V. 
Upper panel: uj = 2.5J; black circle is centered on a pinch- 
point, (b) Section along path segment Vi (see text). 



tributions. 

The stochastic model is microscopic and its main in- 
gredients are a conservation law and the pyrochlore lat- 
tice structure. A long wavelength description is provided 
by the mapping to flux fields 0, [§| and it is interesting 
to see how the dynamics translates under this mapping. 
Taking the low temperature, small q limit, the correla- 
tors for the continuum flux fields B(q, t) implied by the 
stochastic model are 



(i?,(q,i)i?,(-q,0))(x (<5y -*f 



(Ml 



This result can be derived from a Langevin equation for 
the continuum flux fields in which the 'monopole density' 
p = V • B obeys a continuity equation 9fp + V • j = 0, 
with monopole current density 



j = SrATB - 8r Ja^Vp + r^{t) 



(8) 



Here, the second term is the usual diffusion current aris- 
ing from a density gradient, while the first describes re- 
sponse to an entropic force. This response involves a 
drift current of the magnetic charge density p in the field 
B that mimics electrical conduction in electrodynamics 
and is responsible for the flat relaxation rate. Related 
results have been obtained recently in a study of dynam- 
ics in spin ice, represented by the Ising antiferromagnet 
[l5| . In this case monopoles are discrete and it has been 
argued that purely diffusive dynamics are insufficient to 
explain observations and a full description must include 
the network of Dirac strings between monopoles, which 
are essentially entropic, as well as dipolar interactions 
[lit . In spin ice, dipolar interactions lead to Coulomb- 
law forces between monopoles. By contrast, in Eq. [5] 
Coulombic forces appear purely entropically. 



In summary, we have considered wavevector and fre- 
quency resolved dynamics of the classical pyrochlore anti- 
ferromagnet. The relaxational behavior is well captured 
by a stochastic model that conserves total spin. Spin 
diffuses with a diffusion constant independent of temper- 
ature, and entropic forces drive currents to relax con- 
figurations with a rate independent of wave vector and 
inversely proportional to temperature. 
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FIG. 4: (Colour online) Solid blue lines: normalized corre- 
lation function 5(q, t)/S'(q, 0). Dashed red lines: prediction 
from Eq. (O with F = 0.167, both shown at 3 times, t. 



